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Raman and ARPES experiments have demonstrated that in superconducting underdoped cuprates 
nodal and antinodal regions are characterized by two energy scales instead of the one expected in 
BCS. The nodal scale decreases with underdoping while the antinodal one increases. Contrary to 
the behavior expected for an increasing energy scale, the antinodal Raman intensity decreases with 
decreasing doping. Using the Yang, Rice and Zhang (YRZ) model, we show that these features are a 
consequence of the non- conventional nature of the superconducting state in which superconductivity 
and pseudogap correlations are both present and compete for the phase space. 



The pseudogap state of underdoped cuprates is charac- 
terized by a nodal-antinodal dichotomy with Fermi arcs 
at the diagonals of the Brillouin zone (nodal region) and 
a gapped antinodal region 1,2 . Raman^ and angle re- 
solved photoemission (ARPES) experiments 5 have con- 
firmed that a nodal-antinodal dichotomy is also present 
in the superconducting state 6 . These results suggest that 
the superconducting state cannot be simply described by 
BCS theory, contrary to general believe 7 . 

Inelastic Raman scattering probes the zero-momentum 
charge excitations. The response of nodal (xb 29 ) an d 
antinodal regions (xB lg ) can be separated. In the super- 
conducting state pair-breaking peaks appear in the spec- 
tra. As the normal state of cuprates is characterized by 
a not yet understood broad continuum, these peaks are 
better identified in the subtracted response in supercon- 
ducting and normal states Ax Blg , 2g = XbZ, 29 XB lg , 2g - 
With a standard d-wave BCS gap A^ cos(2(/>), as gen- 
erally used for cuprates, the frequency and intensity of 
these peaks in both B2 g and B\ g channels would be only 
controlled by A 5. 

The experiments 3 reveal that in underdoped cuprates 
^XB lg and Axb 29 show pair breaking peaks with op- 
posite evolution with doping, instead of the single en- 
ergy scale. The B\ g peak shifts to higher energy and 
loses intensity with underdoping, while the £>2# peak 
shifts to lower frequency without too much change in 
intensity. The Raman spectrum, specially the intensity 
of the Big peak has been one of the experimental re- 
sults more difficult to understand. A modified BCS gap 
with higher harmonics 3 , and vertex corrections 8 or a very 
strongly anisotropic renormalization of the quasiparticle 3 
have been invoked to explain the B\ g peak frequency and 
intensity, respectively. 

ARPES measurements 5 have also uncovered that un- 
derdoping leads to an increase in the gap in the antinodal 
region A max and a decrease in the slope of the gap at the 
nodes, ^a, resulting in a U-shape dependence of the gap 5 - 
instead of the V-shape expected from d-wave BCS. In the 
standard d-wave model, As gives both and A max . A 
single-energy scale and a linear (V-shape) dependence on 
cos(2</>) of the gap are observed in the ARPES and Ra- 
man spectrum of overdoped cuprates^. 
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FIG. 1: (Color online) Comparison of pseudogap Ar and su- 
perconducting As scales with ARPES nodal va and antinodal 
Amax ones and with the frequency at which B\ g (u>B lg ) and 
B2g (uB 2g ) Raman responses peak. The parameters used de- 
pend on doping as, A R {x)/2 = 0.3(1 - x/0.2), A s (x)/2 = 
0.07(1 - 82.6(x - 0.2) 2 ), t(x) = g t (x) + 0.169/(1 + x) 2 , 
\t'{x)\ = gt(x)\t' \ with t' = -0.3, t"(x) = g t (x)t'o with 
to = 0.2. All energies are in units of the bare hopping to- 



In this letter, we show that the appearance of nodal 
and antinodal energy scales and the suppression of in- 
tensity of AxB lg with underdoping are a consequence of 
the non-conventional nature of the superconducting (SC) 
state in which superconductivity and pseudogap (PG) 
correlations are both present and compete for the phase 
space and not just a modified BCS gap. We assume that 
the PG strength, given by A# vanishes at a topological 
quantum critical point (QCP), and that the SC order pa- 
rameter As follows the critical temperature. Fig. 1 and 
the evolution of the Raman intensity with doping x are 
our main results. Only the nodal energy scales follow 
the non-monotonic dependence of the SC order parame- 
ter, being the slope of the gap at the nodes va a good 
measure of A 5. The antinodal energy scale, i.e. the lo- 
cation of the Big Raman peak oob 19 and the maximum 
gap A max in ARPES, is intimately connected with the 
PG. The Raman spectra is calculated at the bubble level. 
The decrease of intensity in the B\ g channel with under- 
doping is not due to vertex corrections, but to the com- 
petition of PG and superconductivity in the antinodal 
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region. The renormalization of the quasiparticle weight 
g t just enhances this effect. 

We use the ansatz proposed by Yang, Rice and Zhang 9 
(YRZ) for the PG. The YRZ model assumes that the PG 
can be described as a doped spin liquid and proposes a 
phenomenological Green's function to characterize it. 
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Here £k = eok — &t'{x) cos k x cos k y — 2t // (x)(cos2k x + 
cos2k y ) — n p , eok = — 2t(x)(cos k x + cosk y ) and /i p is 
determined from the Luttinger sum rule. £#(k, a;) = 
A^(k) 2 /(cc; + eok) diverges at zero frequency at the umk- 
lapp surface eok, and A#(k) = A j r(x)/2(cos k x — cosk y ). 
The coherent part is similar to the BCS diagonal Green's 
function with the non- trivial difference that in BCS, the 
self-energy diverges at the Fermi surface (FS) and not at 
the umklapp one. Besides there is no off-diagonal com- 
ponent of the Green's function in our case and Ar does 
not break any symmetry. For finite An the quasipar- 
ticle peak is split into two and the FS consists of hole 
pockets close to (±7r/2, ±7r/2). At x Cl Ar(x) vanishes at 
a topological transition and a complete FS is recovered. 
Following predictions of mean field theory^ the coherent 
spectral weight factor, g t = 2x/(l + x), decreases with 
underdoping and vanishes at half filling. We use the same 
parameters proposed in the original paper 9 , see Fig. 1. 

Superconductivity is introduced in the standard way 11 
as in Ref. 9 . The diagonal Green's function becomes 

Cgg*(k,o,) = ,. g « (2) 
u - £k - Eij(k,a;) - £s(k,tj) 

E 5 (k,o/) = |A|(k)|/(o; + £(k) + X R (k,-u)) is the SC 
self energy with As(k) = As(x)/2(cosk x — cosk y ) with 
doping dependence given in Fig. 1. Each quasiparticle 
peak splits into four with energies ±E±. 
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The spectral functions A(k, uo) = — 2ImG(k, a;) and 
B(\h,lj) = — 2ImF(k, uo) with F(k, a;) the anomalous 
Green's function are 

A(k )W ) = 5t 7r{(^) 2 ^ + J E; 1 7) + (%) 2 ^-^k) 

+( w +) 2 ,5( w + i?+) + ( w +) 2 ( 5( w -£+)}, 

B(k,w) = <?fr{WkVk(^ + -Ek) + l5 ( w -- E k)) 

+u+v+(J(o; + ££)+<J(w -E+))}, (3) 

with coherence factors v^ 1 = \ — 6 k / E^ ^ and = 

3 (*k whereat = |(l±(^-^ +A| k )/^) 

and 6± = f k a£ ± A^ k (^ k - fko). 

In the bubble approximation 12 the Raman response is 

lm{ X ,Am = E(^) 2 / J (»fM - + fi)) 

k ^ 

{A(k, + fi)A(k, cj) - £(k, + Sl)B(k, uj)}. (4) 



FIG. 2: Insets in (a) and (b): Raman vertices (7 k l9 ) 2 and 
(7 k 29 ) 2 in the first quadrant of the Brillouin zone (BZ), with 
the Fermi pocket for x = 0.14 drawn. Inset in (c): bands 
for x = 0.14 in the SC state. The arrows signal the possi- 
ble optical transitions. Main figures (a) and (b) subtracted 
Raman response in B\ g and B^ g polarizations divided by the 
coherent weight factor g^, for x = 0.14 and x = 0.18 in the 
underdoped regime and to x c = 0.20. (c) and (d) full Ra- 
man response in the B\ g and B^ g channels respectively for 
the same dopings. Raman units are the same for all dopings. 
The arrows indicate the position, coB lg and cuB 2g i and inten- 
sity of the pair breaking peaks. The 5 functions in Eq.[3]have 
been replaced by Lorenzians with a width of 0.05£(x) with 
t(x) given in Fig. 2a. 



Here uf(oj) is the Fermi function, v = B\ gi B2 g and 7 k , 
the Big and B2 g Raman vertices, are proportional to 
cos k x — cos k y and sin k x sin k y respectively. 

The dependence on x of the nodal and antinodal sub- 
tracted Raman spectra A% can be seen in Fig. 2a and 
2b. Two different energy scales ujb 19 and ojb 29 appear 
for finite Ar. These scales, signalled by an arrow in 
Fig. 2, are plotted in Fig. 1. At x c = 0.2 the PG van- 
ishes and Axbi 9 and Axb 29 peaks appear at an energy 
_close to 2 As-. This is the expected behavior for a BCS 
)superconductor 13 . On the contrary, as x is reduced the 
A\B lg peak shifts to larger frequency and its intensity 
decreases, while the AxB 2g peak shifts to lower frequency 
with a weakly ^-dependent intensity. This behavior re- 
sembles very much the experimental one 3 -. In Fig. 2a 
and 2b the spectra have been divided by g 2 to emphasize 
that the suppression of the intensity in the B\ g chan- 
nel with underdoping is not only due to the reduction 
of the coherent part, as proposed in 3 . With g\ included 
the weakening of the B\ g transition with underdoping is 
enhanced and the B2 g signal decreases. 

The appearance of two different energy scales for B\ g 
and B2 g response is associated to the two pair-breaking 
transitions in the inset of Fig. 2c with energies 2E±(k). 
The distinction between nodal and antinodal signal has 
its origin in the coherence factors i^_(k)t7^.(k) which 
weight each transition. Shown in Fig. 3a (Fig. 3b) for 
x = 0.14, ii 2 _(k)v 2 _(k) (^^_(k)i' 2 _(k)), weights more heav- 
ily the nodal (antinodal) region. The B2 g and B\ g spec- 
tra are respectively dominated by the transitions with 
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FIG. 3: (a) to (c) coherence factors, in the first quadrant of 
the BZ, corresponding to the pair breaking transitions with 
energies 2E± (see text). The weighted densities of transitions 
N wei9ht (E±) are plotted in (d) to (f). 2A S is marked with a 
dotted line. The added coherence factors and total weighted 
density of transitions for x c = 0.2 (critical doping) are shown 
in (c) and (e) recovering the d-wave BCS result. 



energy 2E- and 2E+. The maxima in the Raman spec- 
tra in Fig. 2a and 2b arise from the peaks in the densities 
of transitions N™* ht = £ k u 2 ± (k)v 2 ± (k)5(uo - 2E ± (k)). 

In Fig. 3d N™^* 1 * peaks at a frequency smaller than 
2 As- ^B 2g depends not only on As but also on the 
length of the Fermi pockets along (tt, 0)-(0, tt) as its value 
comes from the edges of the Fermi pockets. Due to 
the shrinking of the pockets with increasing A#, ooB 2g 
is shifted to lower frequencies with underdoping, even if 
a doping independent As is used. In Fig. 3e, the en- 
ergy at which N^^ ht peaks is determined by a saddle- 
point close to (7r,0), along (tt, 0) — (tt, tt) and is bigger 
than 2 As. In the UD region, this energy is mainly con- 
trolled by Ar. If x > x c , only one of the pair-breaking 
transitions gives a finite contribution for a given k. A 
complete d-wave superconducting gap is recovered and, 
as shown in Fig. 3c it is tracked by the added coher- 
ence factors u\(k)vl(k) + u 2 _(k)v 2 _(k). The maximum 



in N\ 



weight 



is not anymore a maximum of the total den- 



sity NZ i9ht 



Nl 



-N' 



weight 



, plotted in Fig. 3f, but 



both densities of transitions match perfectly. N^ ot 
is now the meaningful quantity. It peaks close to 2As 
recovering the BCS single-energy behavior 13 . 

A third crossing transition, with energy E-(k) + 
E+(k), and larger intensity in the B\ g channel is also 
allowed if Ar is finite, as shown in the inset of Fig. 2c. 
Its effect is small in A%, as it is expected in both the PG 
and the SC state. The total responses Xb^ 9 anc ^ Xb% 



including the contribution of this crossing transition in 
the SC state are plotted in Figs. 3c and 3d. It is not easy 
to distinguish this transition from the pair-breaking ones. 
To the best of our knowledge this transition has not been 
found yet in the PG state but it could be hidden in the 
broad background due to strong inelastic scattering. 

As expected for a gap with nodes, the E>2 g response 
in Fig. 3d is linear at low frequencies. The slope is 
doping independent. This independence, observed also 
experimentally 3 , comes from a cancellation between the 
dependencies of the quasiparticle weight squared g 2 , the 
SC order parameter, A^, and the density of states, via 
the renormalization of the band parameters. 

The nodal and antinodal scales can be also seen in 
ARPES. As discussed in ref£ for x < x c in the PG state, 
the FS consists of closed hole pockets. However, due to 
the weak spectral weight of the quasiparticle pole in the 
outer edge of the pocket, the ARPES spectra resembles 
the Fermi arcs observed experimentally. The length of 
these arcs increases with doping, as seen in Fig. 4a and 
Fig. 4b. A complete FS is recovered when Ar = in 
Fig. 4c. In the absence of a complete FS, to analyze 
the k-dependence of the gap we take the surface with 
maximal intensity uo = and A^ = 0. This surface, 
signalled in blue in the pictures, resembles the one inter- 
preted experimentally as the underlying FS. To compare 
with experiments^ we define va as the derivative of the 
energy with respect to cos k x — cos k y at the nodes 14 and 
^max as the maximum gap along this surface. Shown 
in Fig. 4d for x = 0.05, when As is zero but Ar fi- 
nite, the energy vanishes along the arc and a gap opens 
linearly with cosk x — cosk y from the arc edge. A max 
increases with underdoping. A finite A^ opens a gap 
along the arc in Fig. 4e. This gap depends linearly on 
cosk x — cos k y , with slope va very close to A^. Outside 
the aroi^, the gap depends on both Ar and A^. In this 
UD SC region va increases with doping and A max de- 
creases (see Fig. 1). Correspondingly, the spectra does 
not depend linearly on cos k x — cos k y but has a U-shape 
dependence on cosk x — cosk y with a kink around the 
arc edge. Deviations from linearity increase with under- 
doping. In Fig. 4f, at x = x c , the linear V-shape BCS 
dependence re-appears and va and A max converge. 



The evolution of ARPES scales va and A^ 



with 



doping is plotted in Fig. 1, and compared with the ones 
found in Raman and the input A^ and Ar. The similar- 
ity with experimental data is striking 3 . With the As(x) 
used, ujB 2g and va are non-monotonic on doping. On the 
other hand, the frequency at which B\ g peaks oob 19 fol- 
lows very closely 2A max and both decrease as Ar does. 
Twice the gap value at (7r,0), sometimes compared in 
the literature with uob 19 , is expected to be a bit larger 
than 2A max plotted here. The nodal and antinodal scales 
merge when the pseudogap correlations disappear. 

In conclusion, we have reproduced the deviations from 
BCS in Raman and ARPES experiments in underdoped 
superconducting cuprates. Nodal and antinodal energy 
scales with opposite doping dependence appear in both 
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FIG. 4: (a) to (c) Map of ARPES intensity in the first quad- 
rant of the BZ, at zero frequency and zero As for (a) x = 0.05, 
(b) x = 0.14 and (c) x c = 0.20. In blue, the maximum inten- 
sity surface, (d) to (f) Energy spectrum corresponding to (d) 
x=0.05, (e) x=0.14 and (f) x=0.20, along the surface marked 
in (a) to (c). As is finite in (e) and (f). The 5 functions 
in A(k, uj), have been replaced by lorenzians of width O.OOlto 
and the spectral function convoluted with a gaussian of width 
0.02t (~ 6 - lOmeV) and a temperature T = 0.001t . 



the slope of the gap at the nodes, ua, as measured by 
ARPES, is a good measure of while the maximum 
value of the gap A max arises from an interplay between 
the pseudogap Ar and As. We emphasize that we have 
not tried to fit the experiments but just taken the values 
proposed in the YRZ paper 9 . For the values used, x c and 
optimal doping coincide. If this is not the case^ 6 , it is x c 
which controls the emergence of anomalous behavior. 

Similar two-scale behavior could appear in other mod- 
els with a QCP 17 . Possible differences could be the de- 
creasing spectral weight with underdoping, important for 
the constancy of the slope in I>2g channel, which a priori 
is not expected in other QCP models. In the YRZ ansatz 
the FS is truncated without breaking of symmetry and a 
topological transition happens at x c , in agreement with 
experiments 16 and Dynamical Mean Field Theory 18 . Our 
results suggest that the pseudogap is not a precursor to 
the superconductivity but has a different origin and per- 
sists in the superconducting state and that the antinodal 
scale depends on both the pseudogap and the supercon- 
ducting order parameter. The smooth convergence of the 
antinodal scale with the superconducting order parame- 
ter and a peak in ^XB lg are hard to understand in models 
with separation in k-space in which antinodal quasipar- 
ticles, responsible of the pseudogap, do not participate in 
super conduct ivity-i^. While not included here, we believe 
that the appearance of two energy scales in the SC state 
is robust enough to survive inelastic scattering. 



spectra. The nodal B<i g response, peaks at a frequency 
00B 2g which qualitatively follows the doping dependence 
of the superconducting order parameter As. On the con- 
trary, the energy of the pair breaking transitions in the 
antinodal region, ojb 19 , decreases monotonically with in- 
creasing doping and its intensity decreases with under- 
doping due to the competition between pseudogap and 
superconducting correlations. Twice the maximum value 
of the ARPES gap A max along the maximum intensity 
surface follows very closely ojb 19 - Within this model 
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